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Abstract 

We discuss electron - positron pair - production by super - intense, short laser pulses off 
the physical vacuum state locally deformed by (stripped) nuclei with large nuclear charges. 
Consequences of non - perturbative vacuum polarisation resulting from such a deformation 
are shortly broached. Production probabilities per pulse are calculated. 
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Vacuum fluctuations entail observable effects : local deformation of the translationally invariant 
vacuum state leads to vacuum polarisation which is amenable to experimentation. The classic example 
for such a situation is the Casimir effect where two parallel plates deform the vacuum state and thus 
give rise to an attractive potential. Much more efficient polarisations of the vacuum are produced by 
point - like charges, e.g. stripped nuclei with large Z. 

One of the effects produced by such local Coulomb deformations can be observed in 

heavy - ion collisions : non - perturbative vacuum polarisation prohibits the intrusion of strongly 

bound electron energy levels into the Dirac sea spectrum and thus precludes positron production via 

the 'charged vacuum' mechanism proposed by Greiner et al. [[l]. The recently observed absense of 

such positron emission processes is in accord with this mechanism predicting repulsive forces which 

prevent electron levels from penetrating into the Dirac sea even beyond 'critical' nuclear Coulomb 

charges. 

In this letter we shall address yet another physical consequence of local vacuum 
deformation : multi - photon production of electron - positron pairs in super - intense laser fields. As 
we shall see multi - photon transitions involving 10 5 — 10 6 photons occur with measurable probability 
even at laser intensities near the reach of experimentation. Super - intense, pulsed laser beams, hence, 
might be used as probes in an investigation of the structure of the physical vacuum which in itself is 
an interesting aspect of such pair production processes. 

Conceptually, the calculation of pair - production probabilities follows clear guide - lines. Let us start 
with a short description of non - perturbative vacuum polarisation. To encode the many - body nature 
of the problem we use a quantum field theoretic approach and write 



where <p(x,t) denotes the electron Fermi - field and the notation for the Dirac operator is chosen as 
in for the moment 



is the potential, seen by the electron, produced by a point - like nucleus of charge Z e (in numerical 
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calculations finite size effects are taken into account as e.g. in M ). We expand 

tp(x,t) =Y / * a (x)e~ leat a a , (3) 

a 

where 

H D ^ a = (-iad + f3m + V ex ^j * a = e a V a , (4) 
and obtain for the ( symmetrized ) vacuum energy 

if E ^ E (fi) 

\ e a > — m e Q < — rn / 

i.e. the energy of the state ( N is a normalisation factor ) 

V D = N H 4*o, (6) 

e Q <— m 

is the Fock vacuum. The polarisation charge density is obtained as 



E 



if E i*^)! 2 - E 



Clearly, 

= (8) 
for Vext = 0. Upon switching on the Coulomb potential ( screened at large distances ), we find 

Pvac + 0; (9) 

the total charge, however, is a deformation invariant Q and stays zero modulo integer multiples of e. 
As we pointed out ||, the total charge 

Vvac = J d 3 Xp vac (x) (10) 

jumps by — e at the positron threshold —m and thus yields a repulsive potential preventing the 'about 
to intrude' bound electron level from producing a 'charged vacuum' state by entering the Dirac sea 
continuum (DSC)Q. 

Our next step is now to study the effect of a twofold source of vacuum 

deformation : in addition to the large Coulomb field originating from a nuclear charge we employ a 

1 The real number A introduced in ji| and for mathematical consistency has been set A = (as in all our previous 
numerical calculations). Of course, this expresses the original definition of antiparticles formulated by Dirac. 



super-intense pulsed laser field to induce further vacuum deformation. Introducing now a linearly, in 
z - direction, polarised, pulsed laser field 

t 

dt' A(t') cos(ujt')e z , (11) 

with the pulsed (intensity) 2 A(t) and frequency lo, we have, employing the dipole approximation 
(which can be seen to be valid, comparing the appropriate scales with the laser wavelengths under 
consideration), 

Ze 2 

V ext (x, t) = — — + e X(t)x ■ e z cos(wt). (12) 

\x\ 

The dynamics of pulsed laser interactions is efficiently described in a picture of dressed states which 
now build up a dressed Dirac vacuum. Pulsing the laser beam now induces transitions which we 
propose to observe as e + -e~ pairs. Dressed states can be represented as Floquet states 

(H D -id t )-$ a (x,t) = I a V a (x,t) 

= * a (x,t+ — ) (13) 

which are used to decompose the electron field as in (|J) 

<p(x,t) = J2*<*(x,t)e- i ~ e <* t a a (14) 

a 

thus defining creation and annihilation operators for dressed states. 
The dressed Dirac vacuum state is then given as ^] 

jessed = N JJ ( 15 ) 
e a < — m 

More precisely speaking, the spectrum {e a } is continuous and covers the real axis, bound electron 
states turn into broad resonances whose widths increases with the laser intensity. 
In the dressed state picture, a reaction with a pulsed laser beam is described as a transport process 
in which the state {e Q , is dressed and transported (6| along the pulse envelope A(i); during this 
transport transitions to dressed resonances occur which in the switch - off phase of the pulse are turned 
into occupied bound or scattering electron states : after the pulse we observe a bound or scattered 
electron and a hole in DSC : a positron - electron pair has been produced. Of course, double pair 

2 The states a a ^o are for A — * continuously connected to the eigenstates defined in (W) 



production, leaving two holes and either bound or scattered electrons, takes place, the corresponding 
states will have non - vanishing amplitudes in the pulse - transported state with notably suppressed 
probabilities, of course. 

Two circumstances enhance the probability for pair production to take place and the possibility of 
experimental detections : choosing large Z diminishes the distance of DSC to the lowest electron 
bound state and, hence, favours the multi - photon pair production process; the resonance nature 
of dressed bound electron states entails during the reaction a non - negligible energy band width of 
produced pairs and, thus, an energy - dispersed pair - production spectrum. 

There are methods which aim at the calculation of resonance width and position which cannot be 
expected to lead to reliable results for laser intensities which we envisage. Such a calculation would 
be desirable for a detailed, qualitative discussion of the dynamics of pair-production processes; we 
however have to forgo such an analysis. Nonetheless, we include all these effects implicitly by using 
the equation of motion to follow the time-development of the state vector of the system from the 
beginning to the end of the pulse. 

We calculate the probability per pulse P p for transitions 

laser pulse . rms 
\a >E=-e a >m ► \P >e ff >-rn (lo) 

by numerically solving the initial state problem 

9{x,t) ]t=a0 = \p>E=-e a >m (17) 

for the Dirac equation 

(H D -idt)*(x,t)=0 (18) 

where \a >E=-e a >m is the initial positron scattering state (suitably redefining negative energy state 
quantum numbers Q) at t = just before the switch-on of the laser pulse. The state vector ^f(x,t), 
taken for times t larger than the duration of the laser pulse, is assumed to have, following the discussion 
given above, only single - electron state components with amplitudes a^(t/),e Q > — m besides the 
initial state (0). That is to say that the effective space in which our calculations are performed is 
characterised by the projection operator 

P ef/ = \a >< a\ ]E= _ ca +^|/?>< P\ h> _ m , (19) 



where \a > is the already mentioned positron - scattering state and {\(3 >} are 

single - electron bound and scattering states. A moment's thought shows that this selection of states 
is dictated by the Pauli principle incorporated in the field theoretic picture shortly sketched above. 
Needless to say, by solving in this way the initial value problem in a time interval containing switch-on 
and switch-off of the laser beam one fully describes the dynamics (restricted to one-electron (positron) 
states); in particular the finite width effects discussed above are accounted for. 
Choosing tt as a time shortly after the pulse we finally calculate the probability density 

D _ Xjq, c a =-E X/iAe„>-m \ a a 

"p — 1 ' ^ ' 

l~ia,e a = -E 1 

including the average over (positron) states with energy E. The number n(E) of positrons produced 
in reaction (|l^) to be observed in the interval [E, E + dE] is given as 

n(E) = P p (E)p(E)dE (21) 

where p(E) is the density of states in the DSC. 

To give a quantitative idea of the possibility of detecting electron - positron pairs in experiments we 
present a calculation ^] for a conservative case : we take an ultra short pulse with 20 cycles and a 
photon energy of 5eV and calculate P P {E) for nuclei with Z = 70 and Z — 100. The results are 
shown in Fig.l. Realistic pulses encompass several thousand cycles and the P p (E) , s shown should be 
increased by at least an order of magnitude - we conclude that cross sections, i.e. probability per pulse, 
are within the detection capabilities of experimental setups for sufficiently dense samples of stripped 
high Z nuclei - it is beyond the judgement of the authors to estimate the feasibility of preparing such 
samples. 

The measurement of pair production rates is highly desirable since interesting insight into the structure 
of the physical vacuum could be provided in this way. 



3 Numerical calculations were done using the methods outlined in Q : we used a basis set of maximally 1600 
(j = i,...,^) eigenfunctions of the Coulomb Dirac problem to write the time dependent Dirac equation as set of 
ordinary differential equations which, in turn, were solved by the predictor-corrector routine by Shampine and Gordon 
|^|. Careful checks were employed to ascertain that the dimension of our reference space was sufficiently large (actually 
it was seen that stable results already appeared at a significantly lower dimension) to guarantee numbers of the required 
accuracy. 
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Figure Caption 

Fig. 1 : The probability density P P (E) ( see (|2C|)) as a function of \ ma x, X max being the maximal 
field strength of the pulse X(t), 

irt 

A(t) = X max sm 2 ( — ). (22) 
The parameters in the numerical calculation are 

• T p = 20 (cycles of laser period) 

• lo = 5 eV 

• Z = 70 and Z = 100 

• m — 511 keV 

• E = 513 keV 
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